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ABSTRACT 

We examine the properties of galaxies in the Galaxies and Mass Assembly (GAMA) survey 
located in voids with radii > 10 h~^ Mpc. Utilising the GAMA equatorial survey, 592 void 
galaxies are identified out to z 0.1 brighter than Mr - -18.4, our magnitude completeness 
limit. Using the Wyia vs. [NiiJ/Ha (WHAN) line strength diagnostic diagram, we classify 
their spectra as star forming, AGN, or dominated by old stellar populations. For objects more 
massive than 5x10^ M©, we identify a sample of 26 void galaxies with old stellar populations 
classed as passive and retired galaxies in the WHAN diagnostic diagram, else they lack any 
emission lines in their spectra. When matched to WISE mid-IR photometry, these passive and 
retired galaxies exhibit a range of mid-IR colour, with a number of void galaxies exhibiting 
[4.6] - [12] colours inconsistent with completely quenched stellar populations, with a similar 
spread in colour seen for a randomly drawn non-void comparison sample. We hypothesise that 
a number of these galaxies host obscured star formation, else they are star forming outside 
of their central regions targeted for single fibre spectroscopy. When matched to a randomly 
drawn sample of non-void galaxies, the void and non-void galaxies exhibit similar properties 
in terms of optical and mid-IR colour, morphology, and star formation activity, suggesting 
comparable mass assembly and quenching histories. A trend in mid-IR [4.6] - [12] colour is 
seen, such that both void and non-void galaxies with quenched/passive colours <1.5 typically 
have masses higher than 10^^ M©, where internally driven processes play an increasingly 
important role in galaxy evolution. 
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1 INTRODUCTION 

Redshift surveys reveal a remarkable amount of structure, with 
the majority of galaxies located in groups, clusters, and along 
the filaments linking these massive structures (e.g. Huchra et al. 
1983; York et al. 2000; Colless et al. 2001). However, it is the voids 
that represent most of the volume of the Universe. Spanning tens 
of megaparsecs, these void regions are extremely under-dense, 
with a galaxy density less than 20 per cent of the cosmic mean. 
These voids are not empty, and contain a sizeable galaxy popula¬ 
tion (e.g. Rojas et al. 2004; Croton & Farrar 2008; Pan et al. 2012; 
Kreckel et al. 2012). Nevertheless, A Cold Dark Matter (A-CDM) 
cosmology predicts more dark matter haloes in voids than we cur¬ 
rently observe galaxies (Peebles 2001; Tikhonov & Klypin 2009; 
Peebles & Nusser 2010), and a better understanding the evolution 
of these galaxies will help us better understand this discrepancy. 
Void galaxies of all masses are therefore ideal objects in which to 
examine intrinsic vs. extrinsic effects on galaxy evolution, having 
formed far from the nearest cluster mass dark matter halo. 

Galaxies can be broadly divided into two categories based on 
their optical colours: blue cloud galaxies with ongoing star for¬ 
mation, and red sequence galaxies with colours consistent with a 
non-star forming stellar population, though there is colour overlap 
between these two populations (e.g. Taylor et al. 2015). Clearly at 
some point in their histories all galaxies were star forming- it is 
how and where the quenching of star formation takes place that is 
one of the most active areas of research today. These quenching 
mechanisms fall into two classes: intrinsic (driven by galaxy mass) 
vs. extrinsic (environmentally driven) evolution. 

Galaxies follow a strong morphology-density relation, such 
that early-type galaxies with little or no star formation are pri¬ 
marily found in high density regions of the Universe such as 
groups and clusters, while late-type galaxies dominate in low den¬ 
sity environments such as the field (Dressier 1980; Smith et al. 
2005; Park et al. 2007). Similar environmental trends are seen for 
colour and star formation, such that red galaxies with no star 
formation favour high density environments (Lewis et al. 2002; 
Kauffmann et al. 2004; Ball et al. 2008; Bamford et al. 2009). This 
is reflected in the colour-magnitude relation, such that a strong red- 
sequence is observed for cluster members, with the fraction of red 
galaxies decreasing out to the field. The red cluster galaxies typi¬ 
cally have reduced star formation rates relative to galaxies in the 
field (e.g. Balogh et al. 1997, 1998), along with little to no gas (e.g. 
di Serego Alighieri et al. 2007). 

At a naive first glance, it appears that environment is the 
clear driver in galaxy evolution: when a galaxy enters a high 
density region of the Universe (e.g. a cluster), it is stripped of 
star forming material, and it ceases star formation. Processes 
such as galaxy harassment (Moore et al. 1996), viscous strip¬ 
ping (e.g. Rasmussen et al. 2008; Cluveretal. 2013), strangula¬ 
tion (Larson et al. 1980), tidal stripping, and ram-pressure stripping 
(Gunn & Gott 1972), can all disrupt galaxies and remove their fuel 
for star formation (e.g. Wetzel et al. 2012). However, it has been 
observed that the current level of star formation (or colour) in a 
galaxy is more strongly correlated to galaxy mass than local envi¬ 
ronmental density (e.g. Haines et al. 2007; Wijesinghe et al. 2012, 
Alpaslan et al., sumitted), though conflicting results are also found 
(e.g. Balogh et al. 2004). 

While the above environmental processes certainly play an 
important role in quenching star formation in the satellite galaxy 
population, the situation becomes unclear for central, high mass 
galaxies. Galaxies with high stellar masses are typically red in all 


environments. Indeed, high mass, isolated, early-type galaxies ex¬ 
ist, with examples in the nearby Universe including NGC 3332, 
NGC 5413 and IC 1156 (Colbert et al. 2001). These galaxies are 
optically red, with no evidence of star formation (such as Ha emis¬ 
sion) in their nuclear/central spectra. Mass-quenching through pro¬ 
cesses such as AGN feedback and the virial shock heating of in¬ 
falling cold gas are important at high galaxy mass as these pro¬ 
cesses become more efficient (e.g. Woo et al. 2013). Gas-supply 
can also be heavily depleted during episodes of star formation, trig¬ 
gered during merger episodes such as minor mergers with satellite 
galaxies. 

Separating these extrinsic and intrinsic quenching mecha¬ 
nisms is not trivial. By examining galaxies in voids, we are largely 
excluding the environmental effects present in groups and clus¬ 
ters. Major mergers and interactions between void galaxies are 
expected to be rare (though not entirely absent, see Beygu et al. 
2013 for an example of an interacting system in a void). Like 
other isolated galaxies, void galaxies are expected to build up their 
mass primarily via star formation and minor mergers, resulting in 
discy morphologies. These void galaxies are not entirely cutoff 
from the cosmic web, linked via tendrils joining them to galaxy 
filaments (Sheth & van de Weygaert 2004; Zitrin & Brosch 2008; 
Alpaslan et al. 2014), continually supplying them with the fuel for 
star formation, and increasing the chance of minor mergers. 

Examining the void galaxy population is impossible without 
deep, wide-field redshift surveys. Such surveys are vital not only to 
provide a statistically significant sample of void galaxies, but also 
to identify the void regions themselves via accurate positional and 
distance measurements. In this paper, we examine void galaxies 
in the Galaxy and Mass Assembly (GAMA) survey (Driver et al. 
2009), utilising optical and infrared data. In addition to data at op¬ 
tical wavelengths provided by the GAMA and SDSS surveys, we 
furthermore use mid-IR data to search the GAMA void galaxy pop¬ 
ulation for obscured star formation. Sensitive to the emission from 
dust-reprocessed star formation, data from the mid-IR is crucial in 
establishing if a population of truly passive void galaxies exists (in¬ 
trinsic vs. extrinsic evolution). 

Uncovering the nature of mass quenching first requires ob¬ 
taining a population of well-isolated galaxies with no ongoing star 
formation, spanning the stellar mass range at which the fraction 
of quenched galaxies increases (» 3 x 10*** Mg; Kauffmann et al. 
2003). The most basic way to do this is via the colour-mass re¬ 
lation, selecting galaxies that are passive in terms of their opti¬ 
cal colours (they lie on the red sequence), and that do not ex¬ 
hibit spectral features consistent with ongoing star formation, via 
emission line diagnostics such as the Baldwin, Phillips & Terlevich 
(BPT, Baldwin et al. 1981; Kewley et al. 2006) or Wua vs. [Niij/Ha 
(WHAN, Cid Fernandes et al. 2011) diagrams. 

Examining the stellar populations of these optically passive 
galaxies in other wavelengths can reveal a different picture to the 
optical. For example, galaxies residing on the optical red sequence 
may exhibit a UV excess consistent with star formation in the past 
1 Gyr (Yi et al. 2005; Schawinski et al. 2007; Crossett et al. 2014). 
The mid-IR provides further evidence that not all early-type galax¬ 
ies or optically red galaxies are passively evolving. Clemens et al. 
(2009) show that 32 per cent of early-type galaxies in the Coma 
Cluster are not passive, and their mid-IR colours lie off the tight 
Ks - [16] colour sequence traced by quenched galaxies. Ko et al. 
(2013) match quiescent red-sequence galaxies with no Ho- emission 
to Wide-field Infrared Survey Explorer (WISE) 12 pm photometry, 
and find that 55 per cent of their sample exhibits excess mid-IR 
emission consistent with star formation in the past 2 Gyr. Using 
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a combination of mid-IR and optical colours, we therefore gain a 
better understanding of a galaxy’s stellar population than at optical 
wavelengths alone, and identify truly passive galaxies in both the 
void and general galaxy population. 

The aim of this paper is to examine if the mechanisms respon¬ 
sible for the regulation and truncation of star formation differ as 
a function of large-scale environment. We do this by comparing 
galaxies that reside on the cosmic web in clusters, groups, and fil¬ 
aments, with those that reside in the under-dense voids that sep¬ 
arate these structures (though see Sheth & van de Weygaert 2004 
and Alpaslan et al. 2014 for evidence that void galaxies are, in fact, 
linked to the cosmic web via tendrils). The voids we examine have 
a typical galaxy density 20 per cent of the cosmic mean. We use 
the void galaxy population to search for signs of mass quenching 
in some of the most remote galaxies in the Universe using data from 
GAMA, including GAMA-W/^^ data from Cluver et al. (2014). 
Using a combination of optical (u — r) and mid-IR colours, along 
with emission line diagnostics, we find a population of void galax¬ 
ies with masses > 5 x 10® and both optical and mid-IR colours 
consistent with passively evolving stellar populations. By select¬ 
ing these void galaxies, we are largely removing the environmental 
effects that influence galaxy evolution, and we instead focus on in¬ 
trinsic evolutionary processes. If well-isolated void galaxies with 
high masses are a passively evolving population, they must follow 
the same evolutionary pathways as comparable central galaxies in 
the rest of the Universe (i.e. mass quenching). 

This paper is ordered as follows. The GAMA survey is de¬ 
scribed in Section 2, with our photometry and stellar mass sources 
presented in Section 2.1. Our void galaxy sample, non-void com¬ 
parison sample, and completeness limits are presented in Section 3. 
The colour mass relation for void galaxies is presented in Section 5, 
and the selection of passive galaxies via the WHAN line strength 
diagnostic diagram is presented in Section 5.1. WISE photometry 
is utilised in Section 6 to search for ongoing star formation in opti¬ 
cally passive galaxies, and to search for evidence of mass quench¬ 
ing in the void population in Section 6.4. The properties of the high¬ 
est mass void galaxies (M* > lO'® Mo) are examined in Section 7, 
and we examine the merger histories of isolated galaxies via a com¬ 
parison to the Millennium Simulation in Section 8. We discuss our 
results in Section 9, and conclude in Section 10. 

Throughout this paper we use Ga = 0.7, Q.m = 0.3, and 
H(, = 70 km s“', but for the comparison with Pan et al. (2012) we 
use h in our notation to simplify comparison with the prior litera¬ 
ture. Optical magnitudes are given in the AB system (Oke & Gunn 
1983). The WISE survey is calibrated to the Vega magnitude sys¬ 
tem, and mid-IR photometry is therefore presented in the Vega sys¬ 
tem, allowing for easy comparison to the literature. 


2 THE GALAXY AND MASS ASSEMBLY SURVEY 

The Galaxy and Mass Assembly ii (GAMA ii) survey is a multi¬ 
wavelength photometric and spectroscopic survey that covers three 
equatorial regions centred on o- = 9*', d = -1-0.5° (G09), a = 12*', 
6 = -0.5° (G12), a = 14.5\ 6 = -h 0.5° (G15), along with two 
non-equatorial regions. Each equatorial region covers 12° x 5°. 
The spectroscopic survey targets galaxies to r < 19.8 in the 
G09, G12 and G15 equatorial regions, to a high redshift complete¬ 
ness of > 98 per cent (Liske et al., submitted). The majority of 
spectra in GAMA were taken using the AAOmega spectrograph 
(Saunders et al. 2004) on the 3.9 m Anglo-Australian Telescope at 
the Siding Spring Observatory, with this main spectroscopic sam¬ 


ple supplemented by data from other surveys such as the Sloan 
Digital Sky Survey (SDSS) data release 7 (Abazajian et al. 2009). 
The reduction and analysis of the AAOmega spectra is described in 
Hopkins et al. (2013). 

Further details of the GAMA survey characteristics are given 
in Driver et al. (2011) and Liske et al. (2015), the survey input cat¬ 
alogue is described in Baldry et al. (2010), and the tiling algorithm 
for positioning the AAOmega spectrograph fibres is described in 
Robotham et al. (2010). The reduction and analysis pipeline for the 
AAOmega spectrograph is described in Hopkins et al. (2013), with 
the automated redshift pipeline described in Baldry et al. (2014). 
The optical photometry utilised in this work was obtained from 
SDSS imaging in the u, g, r, i, z bands as described in Hill et al. 
(2011). Additional data products used in this work are the stel¬ 
lar mass catalogue (Taylor etal. 2011), the GAMA environment 
measures catalogue (Brough et al. 2013), the GAMA line strength 
catalogue (Gunawardhana et al. 2013), and the GAMA-VWSE cat¬ 
alogue (Cluver et al. 2014). We also utilise the GAMA Galaxy 
Group Catalogue (Robotham et al. 2011) to remove interloper high 
peculiar velocity objects from the void regions, and examine if mas¬ 
sive void galaxies are the central galaxies in their haloes, else if they 
reside in pairs. We only use sources with reliable redshifts (GAMA 
redshift quality flag nQ > 3, Liske et al. 2015) in this work to en¬ 
sure an accurate determination of environment for all galaxies. 


2.1 Photometry and stellar masses 

Stellar masses for the GAMA ii survey are provided by the cat¬ 
alogue of Taylor et al. (2011) with an update to this catalogue to 
include GAMA ii galaxies. The catalogue furthermore provides ab¬ 
solute magnitudes and colours for GAMA galaxies out to z = 0.65, 
k-corrected to z = 0. The values provided in the stellar mass cat¬ 
alogue are derived via stellar population synthesis modelling of 
the galaxy’s broadband photometry, via comparison to Bruzual 
& Chariot (2003) stellar population synthesis models, assuming a 
Chabrier initial mass function. A Calzetti et al. (2000) dust curve 
is assumed. For full details of these stellar mass calculations, see 
Taylor etal. (2011). 

The values in this catalogue are calculated using aperture pho¬ 
tometry (sEXTRACTOR auto photometry), which may miss a signif¬ 
icant fraction of a galaxy’s light. We therefore apply an aperture 
correction from the GAMA stellar mass catalogue to integrated 
values such as stellar mass and absolute magnitude, to account for 
flux/mass that falls beyond the aperture radius used for SED match¬ 
ing. The GAMA stellar mass catalogue provides the linear ratio 
between each object’s r-band aperture flux and the total r-band 
flux from a Sersic profile fitted out to 10 Re (Tayloretal. 2011; 
Kelvin et al. 2012). 

This flux ratio allows us to clean our samples of galaxies with 
spurious masses which are frequently blended with neighbouring 
galaxies, else their photometry and colours are strongly affected by 
nearby bright stars. We select galaxies where the magnitude cor¬ 
rection is -0.2 < magcorr < 0.1 to ensure we only include galax¬ 
ies with reliably measured stellar masses and colours in our void 
and non-void comparison samples. This unequal cut is applied as 
a galaxy with less flux in its Sersic profile out to 10 Re than in 
the r-band aperture has unreliable photometry. The cuts remove 
11 per cent of the void galaxies prior to magnitude and mass cuts. 
These cuts are essential in defining a sample of galaxies with opti¬ 
cally passive colours consistent with negligible star formation. 
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3 VOID GALAXY IDENTIFICATION 

Cosmological voids are the largest (empty) structures in the Uni¬ 
verse, with radii > 10 /?“* Mpc (Pan et al. 2012). The identification 
of these voids requires large volume surveys with spectroscopic 
distances determined for its target galaxies. Narrow field spectro¬ 
scopic surveys do not have the volume to completely enclose the 
large cosmological voids we are interested in for this work. Due to 
the narrow-field nature of the equatorial GAMA survey, with each 
of the three footprints spanning 12° x 5°, a circular void with a ra¬ 
dius R> 10 Mpc will not be completely contained within a sin¬ 
gle GAMA footprint until z » 0.06. Voids in GAMA may therefore 
have centres that reside outside of the GAMA survey footprints, 
and we would miss them and their galaxies if we select voids us¬ 
ing GAMA only. As targets in the GAMA equatorial regions are 
selected from the Sloan Digital Sky Survey (SDSS) (York et al. 
2000), we are able to use the large continuous survey area of SDSS 
(> 7500 deg^) to trace large scale structure in the nearby Universe, 
and therefore voids which lie in the GAMA survey. As a result, our 
search for void galaxies is confined to the GAMA equatorial (G09, 
G12, and G15) regions. 

The void catalogue of Pan et al. (2012) provides the centres 
and sizes of voids located throughout the SDSS survey region, in¬ 
cluding the GAMA equatorial regions out to z = 0.1. The catalogue 
used heliocentric redshifts. The Pan et al. (2012) catalogue uses the 
Void Finder algorithm (Hoyle & Vogeley 2002) to separate SDSS 
into wall and field galaxies, with their field galaxies defined to be all 
galaxies with a third nearest neighbour distance d > 6.3 Mpc. 
This third nearest neighbour distance selects galaxies located in en¬ 
vironments with < 20 per cent of the mean cosmic density, with 
densities <10 per cent expected in the void centres. The wall galax¬ 
ies are used to trace large scale structure (filaments, groups, and 
clusters), and field galaxies are essentially well isolated. Void struc¬ 
ture is then traced by identifying empty spheres between the wall 
galaxies. The sizes of the void regions are defined by the maximal 
sphere with radius that fills the empty space between the wall 
galaxies. We use the void centres defined in this catalogue to find 
all voids located entirely or partially within the GAMA equatorial 
survey footprints, and then search for all void galaxies out to three- 
quarters of the void radius. Objects with z < 0.002 are excluded to 
avoid Galactic sources contaminating our sample. 

Only galaxies with magnitudes brighter than Mr = -20.09 
were used by Pan et al. (2012) in their separation of wall and field 
galaxies. The sizes of voids always depend on the selection criteria 
used to identify them, and by pushing down the mass/luminosity 
function, voids can be infilled by low mass galaxies. Void size will 
decrease with a fainter magnitude limit due to the weaker spatial 
clustering of low mass relative to high mass galaxies (Norberg et al. 
2002; Zehavi et al. 2011). We illustrate this in Fig. 1, using a 
10 Mpc slice though GAMA to highlight large-scale struc¬ 
ture. The galaxies are split into two luminosity bins: one for galax¬ 
ies with Mr < -20.09 mag, and the other for galaxies with 
-20.09 mag < Mr < -18.4 mag, and their positions in comov¬ 
ing cartesian coordinates plotted. When we extend this luminosity 
range down to M,. = -18.4 (the limiting magnitude at the redshift 
limit z = 0.1 of this work), the size of the voids can be seen to de¬ 
crease as the fainter galaxies are less clustered. We therefore only 
examine void galaxies out to 0.75 R^om to remove the less-clustered 
low mass non-void galaxies from the void galaxy sample. By ex¬ 
cluding these void-edge galaxies, we ensure that we remove the 
effects of large-scale environment from our void galaxy sample. 

To the best of our knowledge, the void regions identified in 


Pan et al. (2012) are bona fide voids with radii > 10 Mpc when 
defined using galaxies brighter than Mr = -20.09 mag. However, 
Alpaslan et al. (2014) identifies fine, low density tendrils that ex¬ 
tend into voids defined by bright galaxies. We assume that all galax¬ 
ies examined here belong to their parent void, though they may be 
connected to regions of higher density by these tendril structures. 

We therefore utilise the GAMA environment measures cat¬ 
alogue (Brough et al. 2013) to examine the local environmental 
density of the void galaxies. To do this, we compare the surface 
density E 5 of the void vs. comparison samples. The surface den¬ 
sity measure is based on the distance to the 5th nearest neighbour 
within a velocity cylinder of ±1000 km s“‘. For the void popula¬ 
tion, we find a mean surface density S 5 = 0.09 ± 0.055 Mpc“^, vs 
E 5 = 0.60 ±0.131 Mpc“^ for the non-void comparison sample. We 
confirm that the void galaxies examined here are therefore located 
in extremely under-dense regions of the Universe. 

Due to SDSS DR7 spectroscopic incompleteness in the G09 
region, this region would be under-dense in the catalogue used to 
create the list of voids in SDSS. We therefore exclude voids and 
galaxies in the G09 region with Deck < 0.0°. This incompleteness 
could result in the identification of artificial voids in regions that 
in fact contain galaxy groups and clusters. We also check our void 
galaxy sample for the effects of the survey edges. There is no pref¬ 
erence for red sequence void galaxies to be located at the edges of 
the survey, i.e. they are evenly distributed. 

3.1 Interloper removal 

Galaxies within groups and clusters can have large peculiar veloc¬ 
ities, making them appear at a higher/lower redshift than the over¬ 
density in which they reside i.e. the “fingers of god” commonly 
observed for cluster galaxies. These peculiar velocities can be suf¬ 
ficient to make a cluster galaxy appear within a void during a sim¬ 
ple radial search around a void centre. Large groups/clusters will 
have the highest spread in the peculiar velocities of their member 
galaxies (e.g. Ruel et al. 2014). 

To remove such galaxies from the void sample, we use an 
update of the GAMA galaxy group catalogue of Robotham et al. 
(2011) to identify group galaxies. However, we cannot just remove 
all void galaxies found to reside in groups from our catalogues. 
Simulations predict dark matter substructure and filaments within 
void regions, consistent with a hierarchical model of galaxy assem¬ 
bly (Tikhonov & Klypin 2009; Kreckel et al. 2011b; Rieder et al. 
2013; Aragon-Calvo & Szalay 2013). Evidence for this substruc¬ 
ture has been found, with a small galaxy group identified in a void 
consisting of three galaxies embedded in a common Hi envelope 
(Beygu et al. 2013), hypothesised to be the assembly of a filament 
in a void. To ensure we do not remove such groups, we set a group 
size limit to separate small void groups from interlopers. All void 
galaxies in groups with >10 members are excluded from our final 
catalogue. This limit is selected to keep groups of 2-3 bright galax¬ 
ies and any satellites. Using this method, 16 galaxies with masses 
>10® Mq are identified that belong to a single cluster in G09 with 
a mass 1.6 x lO'"^ Mq, and a velocity dispersion 558 ± 50 km s“' 
(Robotham et al. 2011). These 16 galaxies are removed from the 
void galaxy sample. 

The remaining galaxies all reside in groups of 6 or fewer 
members, else in pairs or isolation. For the remaining groups with 
more than three members, the majority of group members have 
Mr > -20.09 and are less massive than 5 x 10® Mq. Due to their 
low absolute magnitude, these faint galaxies would be absent in the 
wall/field sample of Pan et al. (2012). Indeed, the sum of the lumi- 
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Figure 1. The effect of limiting magnitude on void size. Galaxies with M,. < -20.09 are more strongly clustered than those with —18.00 < < -20.09, and 

voids therefore appear larger when using a brighter magnitude cut to define their limits. 


nosities of all members in the only void group with six members 
(GAMA GroupID 300360) is = -21.4, fainter than a typical 
BCG. Nevertheless, we retain these galaxies in our sample, as they 
are well isolated from the filaments that trace large scale structure. 

3.2 Non-void comparison galaxies 

To allow the comparison of our void galaxy sample to the general 
galaxy population, we construct a comparison sample of non-void 
galaxies in the same volume as the GAMA void galaxies. Only 
galaxies in the G09, G12 and G15 out to z = 0.1024 (the redshift 
of the most distant void galaxy) were included in this comparison 
sample. Our non-void comparison sample contains galaxies from 
a range of environmental density, ranging from clusters through to 
isolated galaxies in filaments, though void galaxies are removed 
from the comparison sample. These non-void galaxies reside in a 
range of local galaxy density, ranging from a surface density E 5 < 
0.05 Mpc“^ (extremely under-dense), through to over-dense regions 
with S 5 > 100 Mpc“^ (i.e. clusters). 

As voids are not purely spherical, we furthermore need to 
avoid the contamination of the comparison sample with void galax¬ 


ies. A void may extend beyond the edges of the maximal sphere 
used to define it if the void does not have a spherical edge. 
We therefore exclude galaxies to a distance of 1 Mpc from 
edges of the maximal spheres used to define the void radii (i.e. 
RtoM + 1 Mpc). This provides a comparison sample of 14,233 
non-void galaxies in the same volume as the void galaxies prior 
to completeness cuts. The distribution of the non-void comparison 
sample is shown Figure 2 as pale blue dots, with the void galaxies 
shown as purple dots. Also included in Fig. 2 are the positions of 
galaxies in SDSS to highlight the presence of large scale structure 
and voids in the GAMA equatorial regions. 

3.3 Completeness 

While we are interested in how the star forming properties of void 
galaxies vary as a function of galaxy mass, spectroscopic complete¬ 
ness cuts based on galaxy mass alone have a strong luminosity bias 
as blue, star forming galaxies are brighter at optical wavelengths 
than passive galaxies for a given stellar mass. As such, the low mass 
end of the galaxy sample will be biased towards star forming galax¬ 
ies. To avoid this luminosity bias in our sample, our completeness 
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Figure 2. Large-scale structure in the GAMA fields (blue points) overlaid on an 10 /i * Mpc slice of SDSS. The void galaxies identified in this work are shown 
as purple points. It can be cleaidy seen that the void galaxies occupy the under-densities in the large scale structure traced by GAMA and SDSS. 


cuts placed on our sample are made using r-band absolute magni¬ 
tudes, with aperture corrections applied (see Section 2.1), to ensure 
no colour bias in our sample. 

The spectroscopic limit of GAMA ii is r = 19.8, and the red- 
shift completeness at this magnitude limit for the GAMA equato¬ 
rial regions (G09, G12, and G15) is 98.5 per cent integrated over 
all magnitudes, and we therefore use this apparent magnitude as a 
limit. 

Our magnitude completeness cuts are shown in Fig. 3. Both 
the void and non-void galaxy populations are shown, along with the 
r = 19.8 spectroscopic limit curve for GAMA. To z = 0.1024, we 
are complete for galaxies brighter than = -18.4 (the horizontal 
orange line in Fig. 3). We will discuss mass and colour complete¬ 
ness in Sec. 4. 

Low surface brightness galaxies are known to prefer a lower 
density environment to high surface brightness galaxies (e.g. 
Rosenbaum et al. 2009; Galaz et al. 2011). We do not apply com¬ 
pleteness cuts to ensure surface brightness completeness, so we are 
likely incomplete for low surface brightness objects. 

This magnitude limit of Mr = -18.4 excludes dwarf galaxies, 
which have stellar masses < 10* Mg, and are typically fainter than 
Mr = -18. We are unable to use dwarf galaxies in our study of 
mass quenching in void regions. However, Geha et al. (2012) show 
that field dwarf galaxies with no active star formation are extremely 
rare, comprising < 0.06 per cent of galaxies fainter than = -18 
in SDSS Data Release 8. No quenched isolated dwarf galaxies with 
stellar masses < 1.0x10* Mg are found in their study, so including 


these faint galaxies is not vital for a study of mass quenching in low 
density environments. 

We furthermore place an upper limit mass limit of 5 x lO" Mg 
on our void and comparison samples. Galaxies in our samples with 
stellar masses higher than this upper limit are typically found to be 
blended with other sources, or have large background gradients, re¬ 
sulting in their calculated masses being unreliable. While galaxies 
more massive than this exist, we do not expect to find such galaxies 
in voids, as major mergers between the massive galaxies required 
to form them are unlikely in these extremely under-dense regions 
of the Universe. 

In Fig. 3, it can be clearly seen in that void galaxies occupy 
under-densities in the galaxy redshift distribution. This is less clear 
at higher z, as the thickness of the wedge increases, sampling a 
wider range of galaxy density, where voids and groups/clusters may 
overlap in this plot. The number of void galaxies in each GAMA 
equatorial region for the mass limited sample are shown in Table 1. 


4 THE COLOUR-MASS RELATION 

The colour-magnitude/mass diagram is one of the most basic di¬ 
agnostic tools when studying galaxy evolution, used to search for 
trends in galaxy colour with respect to galaxy mass or luminos¬ 
ity. Galaxies follow an approximately bimodal colour distribu¬ 
tion, splitting into the blue cloud and red sequence, dependent on 
whether they are currently forming stars. Here, we use the {u - r) 
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Figure 3. The completeness cut used in this work. The horizontal line is the 
magnitude completeness limit at M,. = -18.4, and the right-hand y-axis is 
the redshift limit of this work, which we take to be the redshift of the most 
distant void galaxy {z = 0.1024). The black curve is the magnitude selection 
limit for GAMA (r = 19.8), which has ~ 98.4 per cent spectroscopic com¬ 
pleteness integrated over all magnitudes. A number of filler objects targeted 
with unused 2dF fibres fall below this selection limit. Void galaxies can be 
seen to occupy the under-densities in the non-void galaxy population. 


Table 1. The distribution of void galaxies with My = -18.4 over the three 
equatorial GAMA fields. Due to SDSS spectroscopic incompleteness in the 
G09 region, we remove void galaxies in this region with d < 0° to remove 
the effect of “false” voids appearing in this under-sampled region. 


Region 

a range 
(J2000.00) 

(5 range 
(J2000.00) 

Void 

Non-void 

G09 

129.0°:141.0° 

0.0°:+3.0° 

144 

1331 

G12 

174.0°:186.0° 

-3.0°:+2.0° 

185 

3861 

G15 

211.5°:223.5° 

-2.0°:-i-3.0° 

263 

2626 


colour-mass relation to compare the void sample and non-void 
comparison sample defined in Sec. 3. 

We present the fc-corrected (u - r) colour-mass diagram for 
the GAMA void galaxies as purple stars in the left-hand panel of 
Fig. 4. The stellar masses and /:-corrected colours are provided by 
Taylor et al. (2011), with an update to include GAMA ii galaxies. 
The rest-frame colours are derived from aperture-matched photom¬ 
etry and SED fits of the galaxies optical colours, and are ^:-corrected 
to z = 0. See Taylor et al. (2011) for further details of the colour 
derivation. Also plotted for comparison as blue dots are GAMA 
non-void comparison galaxies to z » 0.102, the highest redshift 
void galaxy identified in our sample. No mass completeness cuts 
have been applied to the colour-mass diagram, to highlight the fact 
that a large number of void galaxies are low mass, blue objects with 
masses <10 ® Mq: dwarf irregular galaxies that are common in low 
density environments. A red sequence is clearly seen for both the 
void and non-void galaxies, though the majority of void galaxies 
are blue systems with masses < 10'° Mq (i.e. sub-M* galaxies). 


At stellar masses >5x10® Mq, the colour-mass relation is clearly 
bimodal, with galaxy colours consistent with both star forming and 
passive galaxies. We therefore take > 5x10® Mq as the mass thresh¬ 
old when comparing the colour and stellar population properties of 
void and non-void galaxies. 

We also plot colour-mass relation using {u - r)i„, colours cor¬ 
rected for internal dust reddening in the right-hand panel of Fig. 4. 
This correction removes galaxies reddened by internal dust from 
the red sequence, again highlighting the fact that the majority of 
void galaxies are blue in colour. In both panels of Fig. 4, a number 
of high mass (> 10'® Mq), red void galaxies with (ii - r) > 1.9 or 
(u - r)j„, >1.6 hint that at least a fraction of void galaxies have 
ceased star formation, despite residing in the most under-dense re¬ 
gions of the Universe. After the dust reddening is applied, 73 out of 
134 of the red sequence galaxies have red colours (u — r)i„, > 1.6, 
consistent with these galaxies hosting quenched stellar populations. 
We will examine the properties of these red galaxies in greater de¬ 
tail to establish their nature. Given their isolated environment, these 
massive void galaxies will be used to establish how galaxy mass (as 
a proxy for halo mass) is responsible for the cessation of star forma¬ 
tion. We first examine the colours and spectra of the void galaxies 
to establish if any have mission/absorption line features consistent 
with quenched stellar populations. 


5 DEFINING A SAMPLE OF PASSIVE GALAXIES 

5.1 Active vs. passively evolving void galaxies 

The red sequence for the general galaxy population contains a mix 
of star forming, active, and quenched galaxies (e.g. Masters et al. 
2010; Tojeiro et al. 2013; Crossed et al. 2014; Taylor et al. 2015). 
To search for signs of ongoing star formation in the red sequence 
population, we use their GAMA/SDSS spectra to identify emission 
lines consistent with ongoing star formation or nuclear activity, and 
utilise their line strengths to quantify the emission mechanism. 

However, this process is complicated by the fact that the pres¬ 
ence of emission lines in a spectrum does not necessarily indicate 
star formation or nuclear activity. Stasinska et al. (2008) show that 
galaxies identified as LINERs on the Baldwin, Phillips & Terlevich 
(1981, BPT) diagram are not necessarily powered by an active nu¬ 
cleus. Instead, old, hot, low-mass evolved stars provide enough ion¬ 
ising photons to mimic nuclear activity (e.g. Binette et al. 1994; 
Stasinska et al. 2008; Yan & Blanton 2012). Galaxies which appear 
to host active nuclei through their classification on the BPT dia¬ 
gram are actually retired galaxies that have ceased star formation, 
with old, evolved stars providing their ionising radiation. These re¬ 
tired galaxies are actually passive galaxies with spectral features 
that mimic LINER activity. For example, using integral field unit 
spectroscopy, Bremer et al. (2013) show that the extended LINER- 
like emission in NGC 5850 is not confined to its nucleus, but is 
distributed over the galaxy. This emission must therefore be a re¬ 
sult of ionisation from distributed sources (likely post-AGB stars), 
rather than a low-luminosity AGN. 

Nuclear activity can also shape a galaxy’s stellar population. 
AGN activity is important in the regulation of star formation, pre¬ 
venting hot gas cooling to form stars, and AGN identified through 
line strength diagnostics are often found on the red sequence or 
in the green valley. We therefore examine the emission/absorption 
lines of the GAMA void galaxy population to identify those with 
active stellar populations dominated by emission lines (star forma¬ 
tion, AGN, LINERs), vs. those with old stellar populations with 
featureless absorption line dominated spectra. 
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Figure 4. Colour-mass relation for GAMA void galaxies (puiple stars). The left-hand panel shows rest-frame (u - r) observed colours, with the rest-frame 
intrinsic (u — r),„f colours (con'ected for internal dust reddening) shown in the right-hand panel. The colours have been /:-corrected to z = 0 in both plots, 
and corrected for foreground extinction. Also shown are non-void galaxies for comparison (blue dots). The majority of void galaxies are low mass, blue, star 
forming irregular and spiral galaxies. A red sequence is seen for both the void and non-void galaxy populations, which remains after the colours have been 
corrected for internal dust reddening. 


The Wyia vs. [Nii]/Hq' (WHAN) diagram (Cid Fernandes et al. 
2011 ) provides a method of separating these different emission 
mechanisms, including those with spectral lines too weak to be in¬ 
cluded in the BPT diagram. It furthermore allows for low-ionisation 
nuclear emission-line galaxies to be separated into weakly active 
AGN, and retired galaxies that have stopped forming stars, with 
hot, low-mass evolved stars providing their source of ionising radi¬ 
ation. Cid Fernandes et al. (2011) find retired and passive galaxies 
to have near identical stellar populations (indeed, occasionally in¬ 
distinguishable), having formed no new stars in the past 100 Myr 
i.e. they have ceased/retired from star formation. 

To separate active from passive/retired galaxies, we construct 
a WHAN diagram for the GAMA void galaxy population using 
the strengths of the Hor and [Nii]/16584, along with the equiva¬ 
lent width of the Hor emission line feature. For full details of the 
GAMA line strength catalogue, see Gunawardhana et al. (2013). 
An update is provided for the GAMA ii sample. The line widths 
and fluxes were measured using a flat continuum, rather than stellar 
population modelling. As a result, partially filled stellar absorption 
lines are provided as absorption, rather than emission, lines, and 
we discuss such galaxies, along with truly passive galaxies, in Sec¬ 
tion 5.1.1 We therefore correct the Ho- fluxes for stellar absorption 
following Hopkins et al. (2003) and Gunawardhana et al. (2013). A 
constant correction of 2.5 A is applied to the equivalent width of the 
Ha line, and the corrected flux is calculated as follows: 


^ _ EWh„ + EWc 

““ EWh„ 


x/h« 


( 1 ) 


where /h„ is the observed Ha flux, and EWc is the constant correc¬ 
tion factor of 2.5 A added to the Ha equivalent width to correct for 
stellar absorption. A minimum ratio of 3 between the equivalent 
width and the error on the equivalent width is required for either 


the Ha or [Nii]/i6584 line for a galaxy to be included in our line 
strength analysis diagram in Figure 5. In future plots, galaxies with 
this ratio < 3 have smaller symbols than those with a ratio > 3. 


We should note here that when examining the optical spectra 
of void galaxies for signs of activity (either through star formation 
or nuclear activity), we are examining single fibre spectra. This 
limits us to either the central 2.1 arcsec of GAMA galaxies with 
spectra obtained using the AAOmega spectrograph, or the central 
3 arcsec of galaxies with SDSS spectra. Galaxies classified as pas¬ 
sive via a line strength analysis of their nuclear spectra may have 
star forming regions outside of their central bulge. If, for example, 
we are sampling a bulge-dominated spiral galaxy, then we will miss 
star formation in its spiral arms, as we are sampling the old, pas¬ 
sive bulge region of the galaxy with an absence of star formation. 
We will explore this issue in Section 6, where we will use mid-IR 
colours to separate centrally passive from truly passive galaxies. 


The WHAN diagram for GAMA void galaxies more massive 
than 5 x 10® Mq is presented in Fig. 5. Galaxies less massive than 
this are primarily star forming, so we exclude them from our plot 
for simplicity. The diagnostic lines of Cid Fernandes et al. (2011) 
are also plotted. Using this diagnostic diagram, galaxies with nu¬ 
clear spectra consistent with passive or retired stellar populations 
have Ha equivalent widths 0 A < EWho- < 3 A. 25 galaxies meet 
this criteria. A large population (170) of star forming galaxies is 
found. Eight weak AGN, and 29 strong AGN are also found in the 
GAMA void sample, though this sample may also contain galaxies 
with composite AGN and star forming stellar populations. Compa¬ 
rable numbers are found for the non-void comparison galaxies. 
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Figure 5. Left-hand panel'. Wua vs. [Niil/Ha (WHAN) diagram (Cid Fernandes et al. 2011) for GAMA void galaxies. While the majority of void galaxies are 
star forming or host AGN, a number of passive/retired galaxies are seen. Right-hand panel: the WHAN diagram for a mass-matched randomly-drawn non-void 
comparison sample. 


5.1.1 Emission-line free galaxies 

The WHAN diagram (and, indeed, any line strength diagnostic di¬ 
agram), does not include galaxies with Hor in absorption, and the 
GAMA line strength classification does not fit stellar population 
models prior to fitting a galaxy. As such, we miss galaxies with 
partially infilled Hor absorption, and those dominated purely by ab¬ 
sorption lines in the WHAN diagram. 

In our sample, sixteen void galaxies have equivalent widths 
measured for their [Nii]/16584 line, but do not have any measured 
Htt emission. Indeed, a visual inspection of their spectra reveals a 
number to have slight infill of the Ha lines due to emission. Eleven 
have [Nii]/I6584 equivalent widths > 0.5 A and we therefore class 
them as retired galaxies in future plots as a heating source is re¬ 
quired to produce this emission. Any remaining galaxies without 
Ha or [Nii]/16584 in emission are classified as passive in future 
plots, and five galaxies meet these criteria. 

This classification scheme results in 40 galaxies more massive 
than 5 x 10^ Mq without a line strength classification. For such 
galaxies, the emission line fitting code was not able to accurately 
measure either their [Nii]/16584 or Ha lines, and hence their stel¬ 
lar population could not be characterised. These galaxies span the 
complete mass range, but the majority have colours (u - r) > 1.9, 
where the galaxy population is dominated by passive/retired galax¬ 
ies or edge-on, dust-reddened disc galaxies. We plot such objects 
in future plots as grey hexagons. 

5.2 The colour-mass relation for galaxies with WHAN 
classifications 

Do the colours of the void galaxies reflect their line strength mea¬ 
surements? To answer this, in Figure 6 we reconstruct the (u - r) 
colour-mass diagram for void galaxies based on their classifica¬ 
tion on the WHAN diagram. Galaxies with Ha in absorption that 
could not be included in the WHAN diagram are added as pas¬ 


sive/retired galaxies, depending on the strength of their [Nii]/t6583 
line. The top left-hand panel shows their observed (u - r) colours, 
and the (u - r) colours presented in the top right-hand panel have 
been corrected for internal dust reddening. The colours have been 
^-corrected to z = 0 in the plots, and corrected for foreground ex¬ 
tinction. When corrected for intrinsic dust reddening, strong AGN 
are almost completely removed from the red sequence, with only 
red and dead void galaxies having retired or passive stellar popula¬ 
tions remaining. 

As a comparison, we include a sample of randomly drawn, 
matched non-void galaxies. For each void galaxy, we select a non¬ 
void comparison galaxy of similar mass (±20 per cent M*), rest- 
frame, colour uncorrected for dust reddening ({u - r) ± 0.15), and 
redshift (z ± 0.01) to ensure we are selecting galaxies at a compara¬ 
ble stage in their evolution. Few high-mass galaxies with AGN-like 
line ratios are seen in both samples, with the red sequence domi¬ 
nated by passive and retired galaxies at M* > 3 x 10*° Mq. 

The distribution of galaxies by WHAN classification on the 
colour-mass diagram is remarkably similar for the void and com¬ 
parison galaxies. A K-S test reveals that for both the void and 
randomly-drawn mass-matched comparison sample, the colour dis¬ 
tributions of the two populations are identical. The p-values are 
>0.1 that both the void and non-void samples are drawn from the 
same population for all five WHAN classes. When we repeat this 
K-S test using galaxy mass rather than colour, both the void and 
the field samples are drawn from the same population for four of 
the WHAN classes, with the exception of the passive galaxies The 
passive galaxies have p = 0.034 that the void and non-void sam¬ 
ples have the same mass distribution. While we cannot completely 
rule out that passive void and non-void galaxies are drawn from 
the sample mass distribution, it can be seen from Fig. 6 that the 
passive comparison galaxies extend to lower galaxy mass than the 
void population. 
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Figure 6. The (u - r) vs. stellar mass diagram for GAMA void galaxies (top row), and for a matched random sample of non-void comparison galaxies (bottom 
row). The symbols have the same meaning as in Fig. 5. The smaller versions of the symbols have an equivalent width vs. equivalent width error ratio < 3. After 
the correction for internal dust reddening, star forming galaxies and galaxies with strong AGN-like features move off the red sequence, and only passive/retired 
galaxies remain. 


5.3 Colour cuts to define passively evolving galaxies 

Using a combination of WHAN line strength diagnostics and op¬ 
tical (u — r) colours, we show in Fig. 6 that the void and non-void 
red sequence is dominated by galaxies with old stellar populations. 
We can therefore use a combination of optical colours and line- 
strength diagnostics to define a colour cut above which the ma¬ 
jority of galaxies are non star forming. We use a lower mass limit 
M^, > 5 X 10® Mq when selecting passive galaxies- below this limit, 
blue galaxies dominate both the void and non-void populations, 
and there is no evidence that isolated galaxies below this limit are 
quenched (Geha et al. 2012). Red galaxies below this mass are typ¬ 


ically dwarf satellites, which are strongly affected by environment 
and must be excluded in a study of mass quenching. 

We identify a. {u — r) colour cut consistent with a passively 
evolving stellar population from the comparison sample defined in 
Section 3.2. The non-void comparison sample is used due to its 
larger sample size vs. the void galaxy sample. First, we identify 
passive vs. star forming galaxies using line strength diagnostics. All 
comparison galaxies more massive than 5 x lO'" M© are classified 
according to their location on the WHAN line strength diagnos¬ 
tic diagram. This analysis is described in full in Sec. 5.1. Galaxies 
with line strengths that quantify them as AGN are removed from 
the sample. This leaves two remaining classes of galaxy: star form- 
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Figure 7. Histogram showing the fraction of star forming and passive galax¬ 
ies as a function of (u - r) colour for the non-void galaxy population. The 
galaxies were selected via their WHAN diagnostic into star forming or pas¬ 
sive. The passive sample includes retired galaxies with hot, old, low mass 
stars as their ionisation source. We exclude AGN from this plot for simplic¬ 
ity. The fraction of star forming galaxies drops as {ii - r) becomes redder. 

ing galaxies dominated by a young stellar population, and pas¬ 
sive/retired galaxies whose optical spectra are dominated by old 
stars. We then construct a histogram showing the (u - r) colour dis¬ 
tributions of the two samples (Fig. 7). It can be seen from Fig. 7 that 
for colours (u-r) > 1.9, passive/retired galaxies dominate in terms 
of number. This distribution in colour by separating star forming 
vs. passive galaxies is similar to the bimodal colour distribution of 
lower luminosity galaxies in Hoyle et al. (2012). Low luminosity 
galaxies dominate both samples in terms of number. We therefore 
take (u-r) = 1.9 as the lower limit of galaxy colour when defining 
a sample of red void galaxies with optical colours consistent with a 
passively evolving stellar population. 

Selecting passive galaxies based on colour alone is not per¬ 
fect, and it can been seen from Fig. 7 that passive and star form¬ 
ing galaxies overlap in colour when (u-r) > 1.5. This over¬ 
lap of passive and star forming galaxies remains after a model- 
dependent correction for intrinsic dust reddening has been applied, 
and we therefore choose to use uncorrected (ii - r) colours when 
selecting our red comparison and void samples. This separation 
may become more complicated when examining void vs. non-void 
galaxies. Hoyle et al. (2012) show that at a given luminosity, void 
galaxies are typically bluer than their non-void counterparts by 
(u-r) » 0.1. However, our cut ensures that we will primarily be 
comparing galaxies that are non-star forming in both the void and 
non-void samples. 

5.4 Obscured star formation 

Searching for star formation at optical wavelengths has drawbacks. 
Dust obscuration can mask low levels of star formation, such that 
the object’s spectrum will appear passive. To examine if obscured 
star formation is present in the void red-sequence galaxy popula¬ 
tion, we therefore go on to examine their mid-IR properties. The 
mid-IR is a more sensitive tracer of recent star formation than the 
optical. PAH emission is a typical feature of star forming galaxies. 
Longer than 8 yum, emission from dust heated by younger stars be¬ 


gins to trace star formation. The optical colours of the galaxies will 
redden and they will join the red sequence within 1-2 Gyr of the 
cessation of star formation. The heating of circumstellar envelopes 
is the main source of mid-IR emission in galaxies, and this is sen¬ 
sitive to star formation over relatively long timescales (> 1 Gyr). 

The optical red sequence contains not only genuinely passive, 
non-star forming galaxies, but a number of late-type galaxies with 
their colours reddened by dust extinction or with a low level of star 
formation which is not sufficient to move them to the blue cloud. 
This low level of star formation would not be picked up in the fibre 
spectroscopy we examine in this work. We instead use WISE mid- 
IR colours as a tracer of recent star formation in galaxies classed 
as passive/retired based on colour and line-strength diagnostics in 
Section 5.3. 


6 GAMA-WISE 

The Wide-field Infrared Survey Explorer (WISE) telescope provides 
this mid-IR data for the GAMA survey. The 3.4 pm (WI) and 
4.6 pm (W2) WISE bands trace the continuum light from evolved 
stars. Wl is most sensitive to stellar light, and W2 is also sensi¬ 
tive to hot dust. W1-W2 is therefore a good colour for identify¬ 
ing galaxies dominated by AGN emission (Jarrett et al. 2011). The 
12 pm W3 band traces the 9.7 ^m silicate absorption feature, as 
well as 11.3 pm PAH and Nell emission line. The W4 band traces 
the warm dust continuum at 22 pm, and is used to trace AGN ac¬ 
tivity and reprocessed radiation from star formation. Thus by ex¬ 
amining the colours of galaxies in the mid-IR, we can compare the 
recent star formation histories of void vs. non-void galaxies (e.g. 
Jarrett et al. 2013). 

In particular, we are interested in revealing star formation in 
red galaxies via the W3 flux (12.0 pm) which will highlight void 
galaxies with current star formation. Donoso et al. (2012) find that 
80 per cent of the 12 pm emission in star forming galaxies is pro¬ 
duced by stellar populations younger than 0.6 Gyr, WISE is there¬ 
fore ideal for identifying galaxies with low levels of nuclear activity 
and star formation that are not easily found in their optical spectra. 


6.1 WISE photometry 

GAMA sources with WISE photometry are identified by cross¬ 
matching the WISE All-Sky Catalogue to the GAMA II observed 
sources catalogue using a 3 arcsec cone search radius. 86 per cent 
of GAMA sources in G09 were detected in the WISE All-Sky sur¬ 
vey, with 82 per cent of G12 sources and 89 per cent of G15 sources 
detected. The data product and its reduction is described in full de¬ 
tail in Cluver et al. (2014). Here, we match the optically selected 
void galaxy sample defined in Section 3 to the GAMA-U75F cat¬ 
alogue using the galaxies GAMA catalogue IDs. 527 out of 577 
void galaxies with masses > 10^ Mg are matched to the GAMA- 
WISE photometry catalogue. Eor the comparison sample, 7233 out 
of 7718 matches are found. 

The majority (70 per cent) of non-matches in both the void 
and comparison samples have masses < 5 x 10® Mg- 10 per cent of 
all galaxies below this mass limit. We therefore exclude all galaxies 
less massive than this mass due to incompleteness, leaving 295 of 
the R75£-void galaxy matches and 4978 non-void galaxies, more 
massive than 5 x 10® Mg in our WISE sample. 

For the void sample, 10 non-matches with masses > 5 x 
10® Mg vs. 295 matched galaxies are found, approximately 
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3 per cent of the sample above this mass limit. All are less mas¬ 
sive than 1.6X 10'® Mq and have (u-r) > 1.6. The non-matches in 
the comparison sample with masses > 5 x 10® Mq exhibit a wider 
range in mass and colour than the void non-matches. Similar to the 
void galaxy sample, these non-matches are a very small fraction of 
the comparison galaxies (< 3 per cent), and we therefore do not 
expect the absence of the non-matches to have a significant impact 
on our results. 

Fewer than 20 per cent of targets in the GAMA equatorial 
fields have S/N> 2 in the W4 band, and we therefore exclude the 
W4 data from our examination of void galaxy colours, instead util¬ 
ising W3 as a tracer of star formation. For both samples, 30 per cent 
have W3 magnitudes (and therefore [4.6] - [12] colours) that are 
flagged as upper limits or null photometry. Where available, upper 
limits on their W3 photometry from the profile-fitted measurements 
from the WISE All-Sky Catalogue are provided. We include these 
upper limits in our examination of the WISE colour-colour diagram 
in Eigures 8 , 9, and 10 as upper limit arrows. 

In this paper, we present WISE colours in Vega magnitudes, 
the native magnitude system of the WISE dataset. The Vega mag¬ 
nitudes can be transformed to the AB-System using the transfor¬ 
mations of Jarrett et al. (2011). However, care must be taken when 
applying these transformations due to uncertainty in the W3 and 
W4 filter response curves (see Wright et al. 2010, and Brown et al. 
2014a for more information). We therefore choose to work in the 
Vega magnitude system in order to compare the mid-IR colours 
of the GAMA galaxies to the literature. The Vega magnitudes are 
then fc-corrected to z = 0 using the spectral energy distributions of 
Brown et al. (2014b). 

6.2 Optically red void galaxies in the mid-IR 

We investigate if the void and non-void comparison galaxies with 
(u-r) colours consistent with quenched stellar populations are truly 
passive via their mid-IR colours. To do this, we utilise the red se¬ 
quence galaxy samples defined in Section 5.3, with stellar masses 
> 5 X 10® Mq and colours (u-r) > 1.9 consistent with a non 
star forming stellar population. The red galaxies identified in Sec¬ 
tion 5.3 are matched to the GAMA WISE photometry catalogue, 
with 95 void and 2125 non-void galaxies with line strength mea¬ 
surements having mid-IR photometry (31 per cent and 50 per cent 
of the WISE void and comparison samples respectively). 

The [3.4] - [4.6] vs. [4.6] - [12.0] colour-colour diagram for 
the GAMA red sequence void and comparison samples are shown 
in Figure 8. The symbols of the void galaxies in Fig. 8 reflect their 
WHAN diagnostic. Also plotted for comparison is the non-void 
sample as blue dots. Most obvious in this plot is how few of the 
optically red void galaxies with (u - r) > 1.9 have [4.6] - [12.0] 
colour < 1.5. Galaxies with [4.6] - [12.0] < 1.5 are typically 
quenched galaxies such as ellipticals with negligible star forma¬ 
tion, though the red sequence void galaxies exhibit a wide range 
of mid-IR colour. The range of [4.6] - [12.0] colours exhibited 
is large for both the void and non-void central populations, with 
0 < [4.6] - [12.0] < 4.4. This spread in colour reflects different 
levels of ongoing star formation activity for the two populations, 
from quenched galaxies and giant ellipticals ([4.6] - [12.0] < 1.5) 
through to starburst galaxies with [4.6] - [12.0] colours > 3.0. Sev¬ 
eral optically red void galaxies with AGN-like spectral features are 
seen with blue [4.6] - [12.0] colours, and these are likely composite 
objects with active nuclei and ongoing star formation. 

Using their [4.6] - [12.0] colours as a proxy for ongoing/recent 
star formation, we compare the distribution of [4.6] - [12.0] colour 



Figure 8. WISE [4.6] - [12] vs. [3.4] - [4.6] ([W2-W3] vs. [W1-W2]) 
colour-colour diagram for void and non-void comparison galaxies with 
masses > 5 X 10® Mq and passive optical colours (u — r)> 1.9. The smaller 
versions of the symbols have an equivalent width vs. equivalent width error 
ratio < 3. The non-void galaxy population occupying the same region of 
the CMD is also plotted for comparison (blue dots). A range of [4.6] — [12] 
colours is seen for both the void and comparison galaxy populations, con¬ 
sistent with a range of recent star formation activity. 

for the red void and non-void populations to search for differences 
between void galaxies and the general galaxy population. For this 
comparison, we only utilise non-void comparison galaxies in the 
same mass range as the void galaxies, with masses 5 x 10® Mq < 
M* < 1.5 X 10" Mq. a two-sided K-S test gives a p-value= 0.42, 
and we therefore reject the hypothesis that the void and non-void 
galaxy samples are clearly drawn from different populations. An 
identical result is found for the [3.4] - [4.6] colours of the void and 
non-void galaxies, with a p= 0.22. Again, that the two populations 
appear to be drawn from the same colour distribution. 

6.3 WHAN-classifled galaxies in the mid-IR 

The mid-IR bands covered by WISE are particularly sensitive to ob¬ 
scured star formation, or enhanced ISM emission (e.g. from nuclear 
activity), and we might therefore be able to use a combination of 
mid-IR colour and WHAN classification to better understand the 
true nature of the stellar populations in these galaxies. Using the 
emission-line classifications presented in Section 5.1, we examine 
the mid-IR colours of the galaxies with the nature of their stellar 
activity classified using the WHAN diagram. The [4.6] - [12] vs. 
[3.4] - [4.6] colour-colour diagram for the GAMA void galaxy sam¬ 
ple is replotted in the left-hand panel of Fig. 9 for galaxies classified 
according to their line strength ratios on the WHAN diagram. 

Of the 261 void galaxies with WHAN classifications identified 
in Section 5.1 with masses > 5 x 10® Mq, 254 are present in the 
GAMA-VWSE catalogue. From these, 30 have [4.6] - [12] colours 
derived from upper limits on their photometry. For the non-void 
comparison sample, 255 have GAMA-W/SE matches, of which 43 
have [4.6] - [12] colours derived from upper limits on their W3 
photometry. 
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As Fig. 9 illustrates, both the void and non-void samples 
show trends in [4.6] - [12.0] colour dependent on their WHAN 
classification. Passive/retired galaxies tend towards bluer mid-IR 
colours, whereas star forming and AGN classified galaxies have 
redder mid-IR colours. It is clear, however, that passive/retired 
galaxies in both the void and comparison sample do not have mid- 
IR colours that reflect their line strength classifications. Galax¬ 
ies with [4.6] - [12] < 1.5 are consistent with no ongoing star 
formation (Jarrett et al. 2011; Cluveretal. 2014), and we utilise 
this cut to separate passive vs. star forming galaxies. The locus 
of elliptical galaxies is located at [4.6] - [12] » 0.5, with spi- 
ral/starburst galaxies having [4.6] - [12] > 2. Fig. 9 shows that the 
passive/retired galaxies for both the void and comparison samples 
have 0 < [4.6] - [12] <4, spanning a wide range of current star for¬ 
mation activity. As we are using nuclear spectra in our emission line 
diagnostics when selecting these passive/retired galaxies, we may 
miss light from non-central star forming regions in these galaxies, 
else their star formation is obscured in the optical. Galaxies classed 
as star forming via WHAN diagnostics in both the void and com¬ 
parison samples have [4.6] - [12] > 1.5, and lie on the region of the 
[3.4] - [4.6] vs [4.6] - [12] WISE colour-colour diagram occupied 
by spiral/star forming galaxies. The WHAN classifications for star 
forming galaxies therefore accurately reflect their current state of 
star formation activity in the mid-IR. 

Galaxies classed as AGN using WHAN line strength diagnos¬ 
tics have [4.6]-[12] > 2.5, and the colours of these galaxies overlap 
with star forming galaxies. Stern et al. (2012) take [3.4] - [4.6] > 
0.8 as the mid-IR colour criterion for luminous. X-ray selected 
AGN. None of these galaxies in either the void or comparison sam¬ 
ples meet this criterion. They do, however, meet the AGN crite¬ 
ria in Jarrett et al. (2011), where AGN can have [3.4] - [4.6] as 
low as 0.5. Meanwhile, weak seyferts can have colours lower than 
this threshold. These galaxies have values of [4.6] - [12] consis¬ 
tent with them being actively star forming or hosting AGN activity, 
but blue [3.4] - [4.6] ~ 0 inconsistent with them hosting power¬ 
ful AGN, though they may contain weak seyferts. The objects are 
low luminosity AGN or LINERs with high levels of star formation 
i.e. they are composite AGN, producing mid-IR SEDs that overlap 
with non-AGN star forming galaxies. 

6.4 Galaxy Mass and the cessation of star formation 

Massive galaxies in all environments are typically red, with 
quenched star formation. By investigating how galaxy colour 
changes as a function of stellar mass (as a proxy for halo mass 
in isolated galaxies), we can examine the role mass plays in the 
quenching of star formation. As optical colours and line-strength 
diagnostics alone are not sufficient to identify truly passive galax¬ 
ies, as we show in Section 6.3, we utilise the GAMA WISE cata¬ 
logue to examine this effect. At mid-IR wavelengths, the polycyclic 
aromatic hydrocarbon feature at 11.3 yum is excited by UV radiation 
from young stars, and WISE [4.6]-[12] colours are a more sensitive 
diagnostic of recent star formation than optical colours. While the 
warm dust continuum at ~ 22 /rm is a better tracer of star formation, 
few void galaxies have sufficient S/N for detections in this band. 
We utilise this colour diagnostic to search for mass quenching in 
the GAMA void galaxy population, using the cut [4.6] - [12] = 1.5 
from Cluver et al. (2014) to separate passive objects from galaxies 
dominated by star formation. See Cluver et al. (2014) for a more 
thorough discussion of using the WISE 12 pm band and [4.6] - [12] 
colour for star formation diagnostics. 

To search for a relation between galaxy stellar mass and star 


formation activity for the void galaxy population, we examine their 
[3.4] - [4.6] vs. [4.6] - [12] colours binned by mass in Eig. 10. 
The galaxies are split into five mass bins, from M* > 10^ Mq to 
M^, < 5 X 10" Mq. The void galaxies are plotted as purple stars, 
with upper limits on the void galaxy colours are plotted as grey 
stars. Also plotted for comparison is the non-void sample as blue 
dots, again binned by mass. We utilise a colour cut [4.6]-[12] < 1.5 
to identify void galaxies with passive mid-IR colours. 

The colour-colour diagram changes with galaxy mass. For the 
lowest mass bin (< 5x10® Mo), both the void and non-void galaxies 
predominantly have IR colours consistent with star formation, with 
[4.6] - [12] > 1.5. It is this region of the WISE colour-colour dia¬ 
gram in which optically selected blue void galaxies with Ha emis¬ 
sion indicative of star formation reside (Figure 9). There is also a 
spread in the [3.4] - [4.6] colours of objects in this low mass bin, 
showing a spread in the level of hot dust emission in these star 
forming galaxies. This spread in [3.4] - [4.6] colours tightens with 
increasing galaxy mass. 

As the mass of the galaxies increases, the [4.6] - [12] colours 
become increasingly dominated by passive colours, with this trend 
seen for both the void and non-void populations. For galaxies 
more massive than lO" Mq, starburst and spiral galaxies with 
[4.6] - [12] > 3.0 are almost entirely absent, and the colours of 
these galaxies instead overlap with those of elliptical galaxies and 
spirals with low levels of star formation. Regardless of environ¬ 
ment, more massive galaxies exhibit lower levels of star formation. 

Twenty-six void galaxies have passive mid-IR colours, and no 
passive void galaxies with reliable photometry exist in our sample 
with masses < 10*° Mq. These passive, high mass void galaxies 
do not exhibit tidal tails or shells in SDSS imaging (they have not 
undergone recent mergers), and do not exhibit edge-on disc mor¬ 
phologies. Furthermore, these massive void galaxies all have lo¬ 
cal galaxy surface density E 5 < 0.35 Mpc“^, which Brough et al. 
(2013) class as being a low-density environment. They have 5th 
nearest neighbour distances > 5 Mpc, and are well isolated in 
SDSS colour imaging, and are therefore unlikely to be cluster in¬ 
terlopers. Galaxies with stellar masses > 10*° Mq typically occupy 
haloes with masses > 10’^ Mq, where various processes such as 
AGN heating and virial shock heating may prevent gas cooling into 
new stars, quenching star formation (e.g. Dekel & Bimboim 2006; 
Croton et al. 2006). 


7 THE HIGHEST MASS VOID GALAXIES 

The twenty-six passive void galaxies identified in Section 6.4 have 
masses > 10*° Mq, where the mass assembly histories of galax¬ 
ies are increasingly dominated by major mergers (Kauffmann et al. 
2003; Baldry et al. 2004). These major mergers between galaxies of 
similar mass (a; 3 : 1) will transform the morphology of a galaxy 
from disc to bulge dominated, removing structures such as spiral 
arms and bars. Alternatively, given their isolation, these isolated 
galaxies may have increased in mass due to star formation, be¬ 
coming quenched above some critical mass threshold, and would 
retain their discy morphologies. We examine SDSS colour imag¬ 
ing for the 48 most massive GAMA void galaxies (regardless of 
their colour) to search for these features, with this high-mass sam¬ 
ple presented in Fig. 11. These massive void galaxies have colours 
1.70 < (m - r) < 2.57, though only two have (u - r) < 1.9. The ma¬ 
jority have optically passive colours using the colour cut defined 
in Section 5.3. The galaxies are sorted from left to right in order of 
increasing mass, with masses 3.6x 10'° Mq < M* < 1.5 x 10** Mq. 
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Figure 9. Left-hand panel'. WISE [3.4] — [4.6] vs. [4.6] — [12] colour-colour diagram for the GAMA void galaxy sample. The colours of the points correspond 
to their classification on the WHAN diagram presented in Fig. 5. The majority of void galaxies with [4.6] - [12] < 2.5 have optical nuclear spectra consistent 
with retired/passive stellar populations. However, a number of these galaxies likely host obscured star formation or are forming stars outside of their nuclei, 
with [4.6] — [12] > 1.5 colours inconsistent with a truly passive stellar population. Right-hand panel: WISE [3.4] - [4.6] vs. [4.6] - [12] colour-colour diagram 
for a mass-matched randomly-drawn sample of non-void galaxies. 


The void galaxies exhibit a range of morphologies, however 
it is clear from Fig. 11 that a large fraction have disc-like struc¬ 
tures, with bar, ring, and spiral arm features common. A number 
of edge-on discs are seen, as well as ring, disc, and bar features in 
face-on galaxies that would be destroyed during a major merger. 
We hypothesise that based on their colours and morphologies, the 
highest mass void galaxies are red spirals and faded discs that have 
extinguished their gas supply, and have not undergone the inter¬ 
actions/mergers required to transform them into elliptical galaxies 
due to their isolation. 

We furthermore compare the morphologies of the void galaxy 
population to a sample of 48 non-void comparison galaxies. For 
each void galaxy, a comparison galaxy was drawn from the non¬ 
void sample defined in Section 3.2. The comparison galaxy was 
selected to have a similar mass (within 20 per cent of the void 
galaxy’s mass), rest-frame colour (u - r) ± 0.15, and redshift 
(z ± 0.01) to the void galaxy. The comparison galaxies are also 
shown in Fig. 11 , underneath their corresponding non-void compar¬ 
ison galaxy. When the void and comparison samples are examined 
by eye, the morphologies of the two populations are similar, with 
the majority of the two samples being disc-dominated. To quantify 
these similarities, we examine the light distributions of the void and 
non-void comparison samples. 

7.1 Light distributions and bulge-disc decompositions 

The bulge-to-disc ratio of a galaxy tells us much about its mass 
assembly history. A pure disc galaxy will have no bulge compo¬ 
nent, having formed its stellar population entirely via star forma¬ 
tion. Likewise, a pure elliptical galaxy will have little/no disc com¬ 
ponent, having formed in a major merger, disrupting its discy com¬ 
ponent. Galaxies with structures intermediate between the two have 
a bulge+disc structure, and will have undergone a mixture of star 


formation, bar-drive evolution and merging, building up their bulge 
component, while maintaining or reassembling their disc. By ex¬ 
amining the bulge/disc ratios of the void vs. non-void galaxies, we 
are able to see which process dominates the two galaxy popula¬ 
tions. 

We expect the the impact of morphological transformation 
from disc-like to bulge-dominated via major mergers to increase in 
significance with increasing galaxy mass, such that bulge-to-total 
mass ratio and Sersic index n will both increase as stellar mass 
increases (Kelvin et al. 2012). The most massive galaxies in the 
Universe are missing from voids (e.g. the brightest cluster galax¬ 
ies and brightest group galaxies), so here we focus on the most 
massive void galaxies available, which only reach a stellar mass of 
1.5 X 10‘' Mq, compared to masses > 5 x 10‘* Mg expected for 
brightest cluster/group galaxies (e.g. Oliva-Altamirano et al. 2014) 

To quantify the morphologies of the high mass GAMA void 
galaxy sample and examine their light distributions, we utilise 
the GAMA Sersic catalogue (Kelvin et al. 2012), and the SDSS 
bulge/disc decomposition catalogue of Simardetal. (2011). The 
bulge/disc catalogue provides three galaxy fitting models: a pure 
Sersic model, an nj = 4 bulge -l- disc model, and a Sersic (free n*) 
bulge - 1 - disc model for 1.12 million galaxies in the SDSS Legacy 
Survey. For the purpose of this paper, we utilise the ni, = 4 bulge -l- 
disc mode to examine bulge-to-total stellar mass ratios. 

We analyse the light distributions of the void and non-void 
comparison galaxies using the GAMA Sersic photometry cata¬ 
logue, with an update to this catalogue for GAMA ii provided by 
Lange et al. (2015). We compare the Sersic indices of the void 
and comparison samples in the r-band for galaxies with M* > 
5 X 10^ Mq. The mass-matched randomly drawn comparison sam¬ 
ple of 304 galaxies was drawn from the non-void sample using the 
same criteria as the Section 5.2. 

Clear trends are found for the Sersic indices of both the void 
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Figure 10. WISE [3.4] — [4.6] vs. [4.6] - [12] colour-colour diagram for void and non-void comparison galaxies (purple stars and blue dots respectively). Void 
galaxies with upper limits on their W3 photometry are shown as grey stars. Five mass bins are shown, as labelled in the top left-hand corner of each panel. 
The vertical green line is at [3.4] — [4.6] = 1.5, the divide between passive vs. star forming galaxies. The comparison and void samples overlap at all masses. 
Passive void galaxies bluer than [4.6] — [12] = 1.5 have masses M* > 10*® Mq, showing that in the most low density regions of the Universe, a fraction of 
massive galaxies have ceased star formation. Unlike their non-void counterparts, massive void galaxies have not undergone extreme environmental quenching, 
and will be star forming unless they reside in sufficiently massive haloes for intrinsic quenching processes to become efficient. 


and comparison sample as a function of mass. Below a stellar mass 
of 3 X 10*® Mq, both samples favour low values of Sersic n, with a 
median Sersic n = 1.27 ± 0.12 for the void population, compared 
to n = 1.39 ± 0.10 for the non-void population. For galaxies more 
massive than 3 x 10*® Mq, we find n - 3.64 ± 0.22 for the void 
galaxy population, and n = 3.68 ± 0.32 for the non-void compar¬ 
ison sample. Within the errors, no difference is found for the light 
distributions of the void and non-void comparison samples. 

We examine the bulge/disc decompositions of our samples us¬ 
ing the r band ni, = 4 bulge + disc model of Simard et al. (201 1). 
212 void galaxies with mass > 5 x 10® Mg are present in the cata¬ 
logue of Simard et al. (2011), and a mass-matched randomly non¬ 
void comparison sample is drawn with the same sample size. Be¬ 
low a stellar mass 3 x 10'® Mq, 3 x 10'® Mq, both the void and 
comparison galaxies have a mean g-band bulge-to-total mass ratio 
mbidge,ghnioiai,g = 0.08. Combined with low Sersic indices ~ 1, in 
this mass regime, both the void and non-void populations are disc- 


dominated, having little/no bulge component, and have likely built 
up their mass via star formation and minor mergers. 

For galaxies with stellar masses > 3 x 10*® Mq, the g-band 
bulge-to-total mass ratios are near identical for the two population, 
with a median value mbuigegIm,otai,g = 0.48±0.05 for the void galax¬ 
ies, vs. 0.54 ± 0.05 for the comparison sample. No obvious differ¬ 
ences are found for the light distributions of the void and compari¬ 
son samples over the entire mass range examined here. 


8 HOW DO VOID GALAXIES ACQUIRE THEIR MASS? 
THE MAJOR-MERGER HISTORIES OE VOID 
GALAXIES FROM SIMULATIONS 

From Fig. 11, many of the highest mass galaxies and their non¬ 
void matched counterparts exhibit discy morphologies, with rings, 
bars, and spiral structures seen in many of the galaxies. They fur¬ 
thermore have masses > 3 x 10*® Mq, the mass range where major 
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Figure 11. Void galaxies and a mass-matched random comparison galaxy selected to have similar properties to the void galaxy in terms of redshift and colour. 
The galaxies are sorted in order of ascending mass. The green circle in each plot is the SDSS single fibre spectroscopy coverage (3 arcsec diameter), and the 
red bar is 5 kpc in length. The galaxies are ordered from from 3.57 x 10*^ M© (top left) to 1.53 x 10^^ M© (bottom right). The galaxies are located between 
z = 0.042 and z = 0.10, and have colours 1.70 < (w — r) < 2.57. It can be seen that a lai'ge fraction of massive void galaxies have disc-like morphologies, with 
spiral arms, bars, and star formation in their outer regions. 
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Figure 11 - continued 
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Figure 12. Left panel: distribution of galaxies that have undergone a major merger at any point during their evolutionary histories (red points), taken from the 
Millennium Simulation. Right panel: the distribution of galaxies that have undergone major mergers in the past 5 Gyr only (purple points). In both plots, all 
galaxies in a 20 Mpc thick slice of the Millennium simulation at z = 0 are shown as grey points. These grey points are common to both panels. The slice 
measures 250 Mpc X 250 /i“* Mpc. The plotted galaxies have stellar mass > 5 X 10^ Mo. Galaxies with major mergers in the past 5 Gyr typically avoid 
the void regions. 


mergers are hypothesised to have an increasingly important role in 
mass assembly. Unless the merger is sufficiently gas rich, major 
merger activity builds the spheroidal/bulge component of a galaxy, 
and destroys disc features, whereas mass assembly via star forma¬ 
tion and minor mergers is less disruptive to the high mass progeni¬ 
tor’s structure. Given they are typically disc dominated, do massive 
void galaxies acquire their stellar mass via star formation or ma¬ 
jor mergers, and does the void merger fraction differ from that of 
the non-void comparison sample? As it is not possible to trace the 
mass assembly history of a galaxy, we must instead utilise results 
from N-body simulations, combined with semi-analytic models, to 
address this question. Rather than presenting observational results, 
the rest of this section uses these theoretical mass-assembly his¬ 
tories (via merger trees) to compare the major merger histories of 
void vs. non-void galaxies. We utilise the Theoretical Astrophysical 
Observatory (Bernyk et al. 2014) to find the time since last major 
merger for galaxies in the Millennium Simulation (Springel et al. 
2005), separating galaxies by local galaxy density. The galaxies 
are modelled using the Semi-Analytic Galaxy Evolution (SAGE) 
galaxy model (Croton et al. 2006). The code only outputs the time 
since the last major merger, though a galaxy may have undergone 
more than one major merger since its formation. 

To visualise how the major merger histories of void vs. non¬ 
void galaxies compare, we take a 250 Mpc x 250 Mpc x 
20 /j“' Mpc thick slice from the original (500 Mpc)^, z = 0.0 
cube of the Millennium Simulation, and examine the distribution 
of galaxies that have undergone a major merger at any point during 
their evolutionary history. To these, we compare the distribution 
of those galaxies that have undergone a major merger in the past 
5 Gyr. The stellar mass range is selected to match that of the GAMA 
completeness limits of our observed galaxy samples. 

We quantify the local environment of the simulated galaxies 
using a 3rd nearest neighbour approach. A (250 /i“’)^ Mpc vol¬ 


ume of the z = 0 Millennium Simulation containing 288354 galax¬ 
ies with stellar masses 5 x 10® Mq to 5 x 10** Mq was utilised to 
find isolated void galaxies vs. non-isolated wall galaxies. The near¬ 
est neighbour search also included a 10 Mpc extension around the 
(250 /r“*)^ Mpc volume, in case the third nearest neighbour lay out¬ 
side the search volume. We use a third nearest neighbour distance 

> 6.3 hr^ Mpc to identify isolated galaxies. 

In Fig. 12, we compare the distributions of all galaxies that 
have undergone major mergers (those having a mass ratio of 0.3 
or higher between the merger progenitors) at some point in their 
history, vs. those that have undergone a major merger in the past 
5 Gyr. Of the galaxies that have undergone a major merger at some 
point in their evolutionary histories, 99.5 per cent have 3rd nearest 
neighbour distances consistent with them lying in a cluster, group, 
or filament. For major merger in the last 5 Gyr (i.e. if most of their 
mass is assembled at recent times), 99.8 per cent are found along 
regions of large scale structure. At all times, major mergers are 
preferably found outside of the void regions, but a small percentage 
of major mergers (< 1 per cent) take place in void regions. 

For the mass range concerned, major mergers are relatively 
rare in both the void and wall galaxy populations. Approximately 
21 per cent of galaxies with stellar masses higher than 5 x 10*** Mq 
have undergone a major merger during their history, and this 
fraction drops to 10 per cent for galaxies with masses between 
5 X 10® Mq and 5 x 10*** Mq. Of the galaxies that have undergone 
major mergers, one void galaxy with a stellar mass > 5 x 10*** Mo 
has undergone a major merger in the past 5 Gyr, while ~ 11 per cent 
of the wall galaxies with comparable mass have undergone major 
mergers in the past 5 Gyr. 

If we assume that major mergers at any past time are the means 
of morphological transformation from late-type discs to early-type 
ellipticals, we may expect ~ 3 of our void galaxies with masses 

> 5 X 10*** Mq to be early-types, and ~ 14 non-void mass-matched 
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comparison galaxies to have early-type morphologies. If we as¬ 
sume major mergers within the last 5 Gyr are the means of mor¬ 
phological transformation (i.e. discs form around remnants for ear¬ 
lier mergers), then we may expect ~ 0 of our void galaxies to be 
early-type and 7 non-void comparison galaxies to be early-types. 

Both major mergers and high mass galaxies are less common 
in isolation than in filaments, clusters, and groups. Star formation 
via gas accretion must be the dominant mass assembly process for 
the highest mass (> 5 x 10*° Mq) isolated galaxies with third near¬ 
est neighbour distances > 6.3 hr^ Mpc, which explains the disc- 
dominated morphologies of the highest mass void galaxies shown 
in Fig. 11. Major mergers are more common for galaxies in higher 
density regions of the Universe, and they are a significant mass as¬ 
sembly process for galaxies with third nearest neighbour distances 
< 1 Mpc and stellar masses > 5 x 10*° Mq i.e. massive galaxies 
in group scale environments or denser. 

From the model and simulation, we conclude that a sim¬ 
ple model of morphological transformation based on the pres¬ 
ence/absence of major mergers in simulations predicts modest dif¬ 
ferences between the void and comparison galaxy samples matched 
in mass. However, we don’t see evidence for this within our sample 
in Fig. 11 nor in the quantitive measures of Sersic indices presented 
in Sec. 7.1. Despite this lack of ellipticals in both our void and 
field observed samples, we cannot rule out the importance of ma¬ 
jor mergers in mass assembly. We expect the majority of elliptical 
galaxies to be found in clusters and massive galaxy groups. Only 
7.5 per cent of galaxies in the low-z GAMA equatorial regions are 
found in groups with ten or more members, so our randomly-drawn 
non-void comparison sample could be lacking in elliptical galaxies. 


9 DISCUSSION 

Using the GAMA ii equatorial survey regions, we have identified a 
population of void galaxies to z » 0.1 brighter than M, = -18.4. 
At the mass range examined here, little to no difference is found 
between the void and non-void galaxy populations in terms of (u-r) 
colour, WISE mid-IR colours, star formation activity, line strength 
ratios, and morphology. Our results are in close agreement with 
those of Alpaslan et al. (submitted), who find that isolated GAMA 
galaxies show a slight trend for higher IR emission, but have similar 
trends in {u - f) colour and morphology to galaxies in filaments. 

By comparing void with non-void comparison galaxies, we 
have found no evidence that quenching is dependent on large-scale 
environment. Quenched galaxies in voids overlap in morphology 
and mass in both samples (Fig. 10 and Fig. 11). This is consistent 
with quenching mechanisms that are a function of halo mass, such 
as virial shock heating, rather than environmental processes. We 
therefore find that galaxy mass is more important than environment 
in quenching star formation, again similar to the results of Alpaslan 
et al. (submitted). We discuss likely methods for mass quenching 
in the void sample, based on their morphologies, line strengths and 
mid-IR colours. 

9.1 Are optically red void galaxies really quenched? 

We classify the void galaxies according to their location on the 
WHAN diagram, to separate galaxies heated by star formation and 
nuclear activity, from those with old stellar populations. Using this 
diagnostic diagram, galaxies with a Hor equivalent width < 3 A 
have passive stellar populations, else are ionised by hot, evolved 
stars, rather than AGN heating. Cid Fernandes et al. (2011) state 


that both passive and retired galaxies have near-identical stellar 
populations, having retired from forming stars over 10 * yr ago. void 
galaxies have stellar populations inconsistent with either star for¬ 
mation or nuclear activity according to their classification on the 
WHAN diagram. 

However, when we examine the location of the passive/retired 
galaxies on the mid-IR colour-colour diagram, it is clear that a num¬ 
ber of them must host optically obscured star formation, else are 
forming stars outside of the central regions sampled by single-fibre 
spectroscopy. A similar result is found using (u-r) > 1.9 as the 
colour cut for quenched vs. star forming stellar populations: the 
red void galaxies show a range of mid-IR [4.6] - [12] colour in¬ 
consistent with them being a completely quenched galaxy popula¬ 
tion, though void galaxies with passive mid-IR colours are found. 
Current Integral Field Unit (IFU) spectroscopic surveys including 
SAMI (Croom et al. 2012) and MaNGA (Bundy et al. 2015) will 
allow us to examine the nature of these galaxies in more detail. 
These IFU surveys provide multiple spectra across the structure of 
a galaxy to > 1 allowing us to characterise the ages of galaxies 
beyond their nuclear regions. 

9.2 Mass quenching in void galaxies 

Do void galaxies undergo mass quenching? To address this ques¬ 
tion, we use the [4.6] - [12] mid-IR colours of the void galaxies as 
a proxy for current star formation activity. Similar to Jarrett et al. 
(2011) and Cluveretal. (2014) we use [4.6] - [12] = 1.5 as the 
division between star forming and passively evolving stellar popu¬ 
lations. In Fig. 10, we bin the galaxies by mass, and examine how 
the distribution of galaxies on the [3.4] - [4.6] vs. [4.6] - [12] 
colour-colour diagram changes with mass. For both the void and 
non-void galaxy populations, the fraction of galaxies with passive 
[4.6] - [12] colour increases with increasing galaxy mass, and pas¬ 
sive void galaxies have stellar masses > 10*° Mq. 

The GAMA galaxy group catalogue (Robotham et al. 2011) 
reveals these intermediate/high mass void galaxies with M* > 
10 *° Mq to be isolated, else they are the dominant/central galaxy in 
their halo (typically a pair or galaxy triplet). The isolated quenched 
void galaxies galaxies are unaffected by environmental processes 
such as strangulation and ram pressure stripping. Unless their gas 
supply has run out, internally driven evolutionary mechanisms must 
be responsible for the quenching of star formation activity in these 
void objects. 

The quenched void galaxies in our sample have masses ~ 
10 *° Mo, and such galaxies typically occupy dark matter haloes 
with Mhaio > 10*^ Mo (e.g. Behroozi etal. 2013). This halo mass 
of 10 *^ Mo is critical in models of feedback processes such as virial 
shock heating and AGN feedback. Croton & Farrar (2008) exam¬ 
ine the halo mass function for both the void and general galaxy 
population in the 2dF Galaxy Redshift Survey. They find that op¬ 
tically defined red sequence galaxies occupy the most massive 
haloes in all environments, and red void galaxies simply occupy the 
most massive dark matter haloes in voids. They find that transition 
from optically blue-to-red dominated dark matter haloes occurs at 
Myir ~ 10 *^ M 0 for all environments, with halo mass the key driver 
in quenching star formation. 

Given their isolation, and therefore lack of environmental 
quenching mechanisms, void galaxies in halo masses < 10 *^ Mq 
will be star forming. For these less massive haloes, gas is accreted 
with temperature 10“*-10^ K, and this gas collapses to form new 
stars. More massive than this critical dark matter halo mass of 
10 *^ Mq, infalling gas is heated by shocks to near the virial tem- 
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perature of 10® K, and stars can no longer form (Dekel & Birnboim 
2006). The highest mass void galaxies we identify with stellar 
masses > 10 *® Mq reside in dark matter haloes with masses suf¬ 
ficient for this process to become efficient, and thus virial shock 
heating may prevent gas cooling and the formation of new stars. 
Processes such as photoheating from the UV background may 
have suppressed star formation in void dwarf galaxies (Hoeft et al. 
2006), but we are unable to identify such galaxies in this work due 
to incompleteness for low mass galaxies. 

Furthermore, feedback by AGNs is hypothesised to become 
efficient in dark matter haloes with masses > 10’^ Mq. This 
AGN feedback preferentially affects the shock heated medium, 
and prevents it from ever cooling to form stars (Dekel & Birnboim 
2006). Croton & Farrar (2008) suggest radio-mode low-luminosity 
AGN heating as a possible source of the shut-down of star forma¬ 
tion for passive void galaxies. The WISE colours of radio-weak 
AGN are indistinguishable from those of passive/retired galaxies 
(Giirkan et al. 2014), and their spectra may not show signs of nu¬ 
clear activity (their stellar component dominates the spectra) We 
are unable to identify such galaxies via the WHAN diagram or mid- 
IR colours (Fig. 5 and Fig. 9). We therefore consider weak radio¬ 
mode AGN heating as another possible quenching mechanism for 
the passive void galaxies. However, with our existing data, we are 
unable to confirm if these low-luminosity AGN exist in our sample, 
though several void galaxies do meet the mid-IR AGN classifica¬ 
tion criteria of Jarrett et al. (2011). 

The majority of void and comparison galaxies in our sam¬ 
ple exhibit discy morphologies in Fig. 11, so we also consider the 
possibility that minor mergers are responsible for the mass accre¬ 
tion, and eventual quenching, of the void galaxies. Minor mergers 
may play an important role in quenching star formation, while still 
preserving discy galaxy morphology (Birnboim et al. 2007). The 
smooth accretion of material over time via minor mergers and gas 
accretion leave the disc of a galaxy intact, while building galaxy 
mass until the critical halo mass of 10 *^ Mq for virial shock heat¬ 
ing is reached. This is followed by bursty star formation and even¬ 
tual quenching. This process can give rise to red and dead late-type 
galaxies in low mass, isolated haloes, with extended shut-down pe¬ 
riods that last 6-7 Gyr. 

The star formation we detect in the mid-IR colour-colour di¬ 
agram in Fig. 9 for optically passive void and non-void compari¬ 
son galaxies could indeed be due to these isolated galaxies shutting 
down gradually over several Gyr. A combination of minor mergers 
and gas accretion is therefore a likely quenching mechanism for 
void galaxies, though a detailed examination of their star formation 
histories is required to confirm this. 

9.3 Mass-build up in void galaxies 

We show in Figure 11 that the most massive void galaxies have 
disc-like morphologies, similar to their non-void counterparts. 
These disc-like structures suggest that gas accretion and minor 
mergers have been responsible for the mass assembly histories of 
these isolated galaxies. Rojas et al. (2004) state that void galax¬ 
ies are bluer than those in richer environments as those in clusters 
and groups have had their gas supply shut-off. If the void galax¬ 
ies are able to sustain their gas supplies over a longer period of 
time than those in high-density regions, the void galaxies will con¬ 
tinue to form stars after those in groups and clusters have been 
quenched (Rojas et al. 2005). Grogin & Geller (2000) suggest that 
luminous galaxies at higher density have less gas and dust available 
for star formation than their void counterparts, due to their earlier 


formation timescale. As a consequence of this difference in gas/dust 
content, luminous galaxies in voids that undergo interactions with 
nearby companions (e.g. minor mergers) are more able to have star 
formation triggered by such an encounter vs. those in groups and 
clusters, and will be more actively forming stars at the present time. 

Evidence for this sustained gas accretion in void galax¬ 
ies is found by Cluveretal. (2010), Kreckel et al. (2011a) and 
Kreckel et al. (2012) who show that many isolated/void galaxies 
have large reservoirs of Hi gas, which favours secularly driven over 
environmentally driven evolution. Unlike their counterparts in re¬ 
gions of higher density, void galaxies are able to retain these large 
gas reservoirs that would otherwise be stripped due to environmen¬ 
tal processes such as mergers. 

Gas accretion as the primary driver of mass growth is further¬ 
more supported by data from simulations, where major mergers are 
rare in isolated galaxies vs. wall galaxies found in filaments, groups 
and voids (Section 8). This comparison shows that 11 per cent of 
the simulated wall galaxies have undergone a major merger in the 
past 5 Gyr, vs. just 2 per cent of the simulated void galaxies. 

While the transformation of discs into spheroids via major 
mergers is well accepted, major mergers may result in disc dom¬ 
inated galaxies provided the merging galaxies are sufficiently gas 
rich (e.g. Hopkins et al. 2009). Thus in a void or indeed any envi¬ 
ronment where environmental processes are weak, galaxies will re¬ 
tain their gas, and more gas-rich than gas-poor major mergers might 
be expected. Major mergers may therefore play a role in mass¬ 
quenching void galaxies, though interactions between similar mass 
galaxies in voids are expected to be rare (though see Beygu et al. 
2013). The models of Hirschmann et al. (2013) show that roughly 
45% of present day isolated galaxies have undergone at least one 
merger, and minor mergers account for 2/3 of these (with similar 
fractions found for galaxies not in isolation). Thus minor mergers 
are more important than major mergers for the mass assembly of 
isolated galaxies. We therefore consider major mergers as a rela¬ 
tively un-important method of mass build-up in void galaxies com¬ 
pared to minor mergers and gas accretion. 


10 CONCLUSIONS 

• We identify void galaxies with masses > 5x10® Mq and iu-r) 
colours consistent with those of a passively evolving stellar popu¬ 
lation. The WHAN line strength diagram reveals the majority of 
these red void galaxies to have ceased forming stars > 100 Myr 
ago. 

• The mid-IR [3.4] - [4.6] vs. [4.6] - [12] colour-colour diagram 
reveals optically selected red-sequence void galaxies to exhibit a 
range of [4.6] - [12] colour, showing a wide spread in their recent 
star formation activity. 

• Regardless of environment, a trend is seen in mid-IR colour 
such that higher mass galaxies are bluer than lower mass galaxies. 
Mass quenching occurs in the void galaxy population, and massive 
galaxies will rarely form stars regardless of their global environ¬ 
ment. 

• All void galaxies with passive mid-IR colours [4.6] - [12] < 
1.5 have stellar masses > 10*® Mq. A similar trend is seen for the 
non-void comparison galaxies. Regardless of global environment, 
central galaxies in sufficiently large haloes have passive stellar pop¬ 
ulations. 

• SDSS imaging reveals similar morphologies for the massive 
void and non-void galaxies at a similar mass, and the two popu¬ 
lations have similar light distributions as indicated by their Sersic 


MNRAS 000, 1-22 (2015) 


GAMA void galaxies 21 


indices. The massive void galaxies are primarily discy in nature, as 
expected for galaxies in a low density environment. 
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